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ABSTRACT 

The parameter space favourable for the resonant excitation of free oscillation modes by 
dynamic tides in close binary components is explored using qualitative considerations 
to estimate the order of magnitude of the tidal force and the frequency range covered by 
the tidally induced oscillations. The investigation is valid for slowly rotating stars with 
masses in the interval between 2 and 20 M©, and an evolutionary stage ranging from 
the beginning to the end of the main sequence. Oscillation modes with eigenfrequencies 
of the order of five times the inverse of the star's dynamical time scale Td yn - i-e. the 
lowest-order p-modes, the /-mode, and the lowest-order g + -modes - are found to be 
outside the favourable parameter space since their resonant excitation requires orbital 
eccentricities that are too high for the binary to stay detached when the components 
pass through the periastron of their relative orbit. Resonances between dynamic tides 
and <? + -modes with frequencies of the order of half of the inverse of the star's dynamical 
time scale on the other hand are found to be favourable for orbital periods up to 
- 200 r dyn, provided that the binary mass ratio q is larger than 1/3 and the orbital 
eccentricity e is larger than ~ 0.25. This favourable range comes down to orbital 
periods of up to 5 — 12 days in the case of 2 — 20 M Q zero- age main-sequence binary 
components, and orbital periods of up to 21 — 70 days in the case of terminal main- 
sequence binary components. 
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1 INTRODUCTION 

In close binaries with non-synchronously rotating compo- 
nents, each star experiences the time-dependent tidal force 
exerted by its companion. The tidal force gives rise to forced 
nonradial oscillations with frequencies that are determined 
by the orbital period and the rotation rates of the component 
stars. The forcing frequencies may, at any time, come close 
to the eigenfrequencies of the stars' free oscillation modes 
because of secular changes in the orbital parameters due to 
orbital and rotational angular momentum losses or because 
of changes in the eigenfrequencies due to stellar evolution. 
Such resonances lead to an enhanced tidal action and re- 
sult in the excitation of the oscillation mode involved in the 
resonance. 

The resonant excitation of free oscillation modes by the 
tidal action of a companion was first suggested in a sem- 
inal paper by Cowling (1941). Since then, various authors 
have approached the problem from both an analytical and 
a numerical point of view with the aim to explain some of 
the evolutionary aspects and the observational characteris- 
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tics of binary stars (e.g. Zahn 1970, 1975; Savonije & Pa- 
paloizou 1983, 1984; Willems, Van Hoolst & Smeyers 2003, 
and references therein). The subject recently received a new 
impetus, first from the discovery of giant planets in close 
orbits around solar-type stars (see Willems et al. 1997, and 
references therein), and subsequently from the increasingly 
accurate large-scale surveys of pulsating stars in close bina- 
ries (e.g. Harmanec et al. 1997, Aerts et al. 1998, Aerts et 
al. 2000) . The potential of tides as an excitation mechanism 
was furthermore found to be particularly promising for the 
excitation of gravity modes in hot B subdwarfs by Fontaine 
et al. (2003). Existing and upcoming space missions such as 
WIRE, MOST, MONS, COROT, and Eddington are likely 
to contribute even further to the rising interest in tidally 
excited oscillation modes. 

The occurrence of resonances and their effects on the 
evolution and the observational properties of a binary de- 
pend on the period and the eccentricity of the orbit, on the 
masses and the radii of the component stars, and on the 
properties of the oscillation mode involved in the resonance. 
The latter in turn depend sensitively on the internal stel- 
lar structure of the tidally distorted star. The multitude of 
these dependencies and the complexity of their interaction 
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with each other make it difficult to form a concise intuitive 
picture on the possibility of exciting nonradial oscillations 
by tides in close binaries. A systematic study unravelling 
the role of the different stellar and orbital parameters would 
therefore provide a valuable addition to the theory of tidally 
induced oscillations as well as to observational campaigns 
dedicated to the search of free or forced oscillations in bi- 
nary stars. 

In this paper, we set a first step in an attempt to under- 
stand the parameter space of stellar and orbital parameters 
for which the circumstances for the resonant excitation of 
free oscillation modes are favourable. The primary aim is 
to clarify the qualitative behaviour of the acting parameters 
and the interaction between them without detailed calcu- 
lations involving the numerical integration of the systems 
of differential equations governing free and forced oscilla- 
tions in close binary components. The results of this study 
are then to serve as a starting point for detailed numerical 
calculations which will be presented in a follow-up investi- 
gation. 

The plan of the paper is as follows. In Section 1, we 
present the basic assumptions and the basic ingredients of 
the theory of forced oscillations in close binaries. Some par- 
ticular attention is paid to the Fourier decomposition of the 
tide-generating potential, which plays a crucial part in the 
investigation. In Section 2, we present a detailed analysis of 
the Fourier coefficients as functions of the orbital eccentric- 
ity. In Section 3, we examine the ratio of the tidal force to 
the gravity at the star's equator. In Section 4, the results 
of the preceding sections are used to constrain the range of 
forcing frequencies that is favourable for the resonant ex- 
citation of free oscillation modes by dynamic tides in close 
binaries. In the final section, we summarise our results and 
present some concluding remarks. 



2 BASIC CONCEPTS AND ASSUMPTIONS 

Consider a close binary system of stars orbiting each other 
under the influence of their mutual gravitational force. Let 
.Porb be the orbital period, a the semi-major axis, and e the 
orbital eccentricity. The first star, with mass Mi and radius 
Ri, is assumed to be a uniformly rotating main-sequence 
star with a rotation axis perpendicular to the orbital plane. 
The magnitude of the rotational angular velocity Q is as- 
sumed to be small in comparison to the inverse of the star's 
dynamical time scale: 



1 



(i) 



where G is the Newtonian constant of gravitation. The sec- 
ond star, with mass M2, is treated as a point mass. 

Furthermore, let r = (r, 8, <j>) be a system of spherical 
coordinates with respect to an orthogonal frame of reference 
that is corotating with the star, and let et W (r, i) be the 
potential giving rise to the tidal force. Here, et is a small 
dimensionless parameter which is proportional to the ratio 
of the tidal force to the gravity at the star's equator and is 
defined as 



et = 



Ui Mi 



(e.g. Tassoul 1987). 

The tide-generating potential can be expanded in 
terms of unnormalised spherical harmonics Y l [ n (9, <f)) and in 
Fourier series in terms of multiples of the companion's mean 
motion n — 2n/P OI b as 
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where ot — k n + m Q is a forcing angular frequency with 
respect to the corotating frame of reference and r is a time 
of periastron passage. The Fourier coefficients e^ m> & are de- 
termined by 
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where Pj m| (x) is an associated Legendre polynomial of the 
first kind, and M and v are respectively the mean and the 
true anomaly of the companion in its relative orbit (e.g. 
Polfliet & Smeyers 1990; Smeyers, Willems & Van Hoolst 
1998). The Fourier coefficients obey the property of symme- 
try Q_ ra _i = c-i t m,k and are different from zero only for 
even values of i -+- \m\. In the particular case of a binary 
with a circular orbit they are different from zero only for 
k = — jti. 

It follows that, for each degree £ of the spherical har- 
monics Y™(9, cj>), the tide-generating potential generates an 
infinite number of partial dynamic tides with forcing angular 
frequencies gt 7^ 0. When one of these forcing frequencies 
is close to the eigenfrequency of a free oscillation mode, the 
tidal action exerted by the companion is enhanced and the 
oscillation mode is resonantly excited. The occurrence of 
resonances between partial dynamic tides and free oscilla- 
tion modes is particularly relevant for the excitation of free 
oscillation modes g + since, for binaries with short orbital 
periods, their eigenfrequencies may be in the range of the 
forcing frequencies induced by the companion. 

Since the terms associated with the third- and fourth- 
degree spherical harmonics in the expansion of the tide- 
generating potential contain an additional small factor Ri/a 
or (Ri/a) 2 with respect to the second-degree terms, the lat- 
ter ones are usually dominant. For the remainder of the 
paper, we therefore restrict ourselves to the terms associ- 
ated with 1 = 2. The only non-zero Fourier coefficients in 
Expansion JHJ of the tide-generating potential are then the 
coefficients C2,-2,fc, C2,o,fc> and C2,2,fc- These coefficients are 
independent of the ratio R\ja so that they are solely de- 
termined by the orbital eccentricity e. Since they also obey 
the symmetry property C2 t - m ,-k = C2, m ,fe, it is furthermore 
sufficient to consider positive values of k only. 

The effects of resonant dynamic tides on the evolution 
and the observational characteristics of a binary depend on 
the strength of the tide and on the properties of the oscil- 
lation mode involved in the resonance (e.g. Smeyers et al. 
1998, Willems et al. 2003). The strength of the tide is deter- 
mined by the mass ratio q = M2 /Mi and the orbital separa- 
tion a via the dimensionless parameter et, and by the orbital 
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eccentricity e via the Fourier coefficients Ct^ m ,h- The prop- 
erties of the osciliation mode contribute to the effects of the 
resonance via the so-called overlap-integral which is propor- 
tional to the ratio of the work done by the tidal force to the 
mode's kinetic energy (for a precise definition see, for exam- 
ple, Press & Teukolsky 1977, Kumar, Ao & Quataert 1995, 
Smeyers et al. 1998). For main-sequence stars, the behaviour 
of the overlap- integral is such that the coupling between the 
oscillation mode and the tidal force generally weakens with 
increasing radial order of the mode. 

In the following sections, we focus on the role of the 
stellar and orbital parameters affecting the strength and the 
frequency of the tides. The role of the internal structure of 
the star and the properties of the oscillation modes will be 
the subject of a forthcoming investigation. 



the dominant contributions generally stem from the terms 
associated with the coefficients C2,o,fc- 

Although the value of k beyond which the contributions 
to the expansion of the tide-generating potential are consid- 
ered to become negligible depends on the desired accuracy of 
the Fourier decomposition, a reasonable estimate for a cut- 
off value can be obtained by observing the rapid decrease of 
the absolute value of the Fourier coefficients from values of 
the order of 10 -2 to values of the order of 10 -4 . If we focus 
on the coefficients associated with m = —2, the maximum 
fc- value that needs to be considered increases from k w 5— 10 
in the case of the orbital eccentricity e = 0.2 to k ~ 25 — 45 
in the case of the orbital eccentricity e = 0.6. In the latter 
case there are thus up to 5 times more forcing frequencies 
induced in the star which may be resonant with the star's 
free oscillation modes g + . 



3 THE ORBITAL ECCENTRICITY 

The orbital eccentricity contributes to the determination of 
the strength of the partial dynamic tides through the Fourier 
coefficients C2, m ,k- For a given eccentricity, the coefficients 
generally decrease with increasing values of k, but the de- 
crease is slower for higher orbital eccentricities. In addition, 
for a given value of k, the coefficients C2, m ,fc generally in- 
crease with increasing values of e. Hence, there is a finite 
number of non-trivially contributing terms in the expansion 
of the tide-generating potential which increases with increas- 
ing orbital eccentricities. 

In order to estimate the value of k beyond which the 
contributions to the expansion of the tide-generating poten- 
tial become negligible as a function of the orbital eccentric- 
ity, the logarithms of the absolute values of the non-zero 
coefficients C2, m ,fc are displayed in the left-hand panels of 
Figs.^Elas a two-dimensional function of k and e. For clar- 
ity, lines for which |c2, m ,fc(e)| « 10 -2 , 10 -3 , 10~ 4 are added 
to the left-hand panels of the figures and one-dimensional 
variations of the coefficients C2, mj fc as a function of k are 
shown in the right-hand panels in the case of the orbital ec- 
centricities e = 0.4, 0.6, and 0.8. The smallest Fourier coeffi- 
cients considered have absolute values of the order of 10~ n . 

As mentioned in the beginning of this section, the gen- 
eral tendency of the coefficients C2, m ,h is to decrease rapidly 
with increasing values of k. The behaviour is slightly more 
complex for the coefficients C2,_2,fc for which, depending on 
the value of the orbital eccentricity, a local minimum and a 
local maximum may occur prior to the rapid decrease with 
increasing values of fc 1 . This is particularly clear in the cross 
sections shown in the right-hand panel of Fig. Q A similar 
behaviour is observed for the coefficients C2,2,fc, albeit some- 
what less pronounced. For a given value of the orbital ec- 
centricity e and sufficiently large values of the Fourier index 
k, the largest contributions to the expansion of the tide- 
generating potential usually stem from the terms associated 
with the coefficients C2,-2,fc, while the contributions of the 
terms associated with the coefficients C2,2,fc are usually neg- 
ligible. For large orbital eccentricities and small values of k, 



1 Note that the local maximum occurs at values of k close to 
those for which kn ss f!p, with flp = n [(1 + e) / (1 — e) 3 ] 
the orbital angular velocity at the periastron of the relative orbit. 



4 THE PARAMETER e T 

The small dimensionless factor et in the definition of the 
tide-generating potential can be used as a first estimate for 
the order of magnitude of the tidal force exerted on a binary 
component by its companion. In particular, the largest term 
in Expansion (J^J of the tide-generating potential, evaluated 
at the star's surface, is of the order of 
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where the maximum is taken over all admissible values of 
m and fc. For orbital eccentricities less than or equal to 0.8 
the maximum takes values between 0.5 and 2.5 [see Eq. £j) 
and Figs. 0£U ■ The eccentricity therefore affects the order 
of magnitude by less than a factor of 2.5. 

Definition JS) of the small parameter et can be rewrit- 
ten by means of Kepler's third law as 



et 
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The parameter is thus proportional to the square of the 
star's dynamic time scale and inversely proportional to the 
square of the orbital period. In addition, it depends on the 
mass ratio 9 by the factor 9/(1 + 5) which is always smaller 
than 1. In the case of a binary consisting of two 10 Mq zero- 
age main-sequence stars with an orbital period of 6 days, the 
parameter et is of the order of 10 -3 . 

For detached binaries, an upper limit for et can be de- 
rived from the requirement that the radius of the star must 
be smaller than the volume equivalent radius of its Roche 
lobe. The latter may be approximated by means of Eggle- 
ton's (1983) fitting formula 



Ru 



O.499 



-2/3 



0.6 9- 2 / 3 + ln(l + 9- 1 /3) : 



(7) 



where 9 = M2/M1. The relative error of the approximation 
is smaller than 2% for < 5 < 00. By means of Definition 
and Eq. 0, one then derives the inequality 



et < 



0.12 



5[0.6 9- 2 / 3 + ln(l + 9- 1 /3 



(8) 



Taking into account the 2% uncertainty in the fitting for- 
mula given by Eq. Q, it follows that < et < 0.13. 
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Figure 1. Left: The logarithm of the absolute value of the coefficients C2,_2,fc(e) as a two-dimensional function of the Fourier index k and 
the orbital eccentricity e. The solid lines indicate the values of k and e for which |ca _ 2 k{ e )\ ~ 10 -2 , 10 — 3 , 10 — 4 . Right: Cross-section 
of the variation of the coefficients C2 —2 fc(e) for the orbital eccentricities e = 0.4, 0.6, 0.8. 
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Figure 2. As Fig.0 but for the coefficients C2 fc( e )- 



Equation Q for the Roche-lobe radius of a binary com- 
ponent is strictly speaking only valid for binaries with cir- 
cular orbits. In the case of a binary with an eccentric orbit, 
the radius of the Roche lobe varies periodically in time, with 
the smallest value being reached when the companion passes 
through the periastron of its relative orbit. An instantaneous 
estimate for the Roche-lobe radius at a time of periastron 



passage is obtained by replacing the semi-major axis a in 
Eq. with the periastron distance r v = o(l — e). The re- 
quirement that the star fits within its Roche-lobe when the 
companion is located in the periastron of its relative orbit 
then yields 



q[0.6q- 2 / 3 +ln(l + q- 1 / 3 )} 3 ' 
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Figure 4. The maximum values of the small parameter ex for 
which a binary with an orbital eccentricity e and a mass ratio q 
is still detached. 



The variation of the right-hand member of this inequality 
as a function of the mass ratio q is illustrated in Fig. 0] 
for orbital eccentricities ranging from 0.0 to 0.6. It follows 
that, for a given value of e, the parameter et takes values 
in the range < et < 0.13(1 — e) 3 . The maximum value 
of et that can be reached by detached binaries decreases 
by approximately one order of magnitude when e increases 
from to 0.5, and by approximately two orders of magnitude 
when e increases from to 0.8. 



Equations JSJ and © can furthermore be combined to 
derive the shortest orbital period for which a binary is still 
detached when the components are located at the perias- 
tron of their relative orbit. The resulting critical orbital pe- 
riods separating detached from semi-detached systems are 
displayed in Figs. and E] for 2 — 20 Mq main-sequence 
stars at the beginning (ZAMS) and at the end (TMS) of 
the main sequence, respectively. The dashed lines represent 
the critical orbital periods in the case of the binary mass ra- 
tio q = 0.1, while the solid lines represent the critical orbital 
periods in the case of the binary mass ratio q = 1. The grey- 
shaded region in between the lines associated with the two 
mass ratios only serves as a visual aid to identify the lines as- 
sociated with the same orbital eccentricity. We furthermore 
note that, for brevity, we here neglect the short phase of 
core contraction that takes place between the exhaustion of 
hydrogen in the core and the start of the Hertzsprung gap, 
and we refer to the stage where hydrogen is exhausted in the 
core as the end of the main sequence. The ZAMS and TMS 
stellar radii appearing in Eq. © are determined by means 
of the analytic approximation formulae derived by Hurley, 
Pols & Tout (2000). 

Since the critical orbital periods separating detached 
from semi-detached systems are proportional to the star's 
dynamical time scale [see Eqs. © and JjJJ], they increase 
with increasing mass of the star. For a given orbital eccen- 
tricity and a given mass ratio, the same dependency yields 
longer orbital periods for stars on the TMS than for stars on 
the ZAMS. In addition, the critical orbital periods increase 
with increasing values of the mass ratio q because the inner 
Lagrangian point L\ moves closer to the star with mass M\ . 
In the next section, we will use the behaviour of the critical 
periods to constrain the range of forcing frequencies induced 
by the tidal force exerted by the companion. 
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Figure 5. The critical orbital periods separating detached from 
semi-detached systems for 2 — 20 Mq ZAMS stars and for binary 
mass ratios q = 0.1 (dashed lines) and q = 1 (solid lines). The 
grey-shaded region in between the lines associated with the two 
different mass ratios connects the lines associated with the same 
orbital eccentricity. 



5 THE FORCING FREQUENCIES 

In order for a free oscillation mode g + to be resonantly ex- 
cited by a dynamic tide, the forcing frequency of the tide 
must be close to the mode's eigenfrequency. Since the forcing 
frequencies depend on the orbital period and the rotational 
angular velocity of the tidally distorted star, the requirement 
that the forcing frequencies be in the range of the eigenfre- 
quencies of the lower-order <; + -modes imposes constraints on 
the parameter space available for the occurrence resonances. 

The eigenfrequencies of the lower-order second-degree 
g + -modes in 2Mq and 20 Mq main-sequence stars at the 
beginning and at the end of the main sequence are shown in 
Fig. [7] in units of the inverse of the star's dynamical time 
scale. The stellar models consist of a convective core and a 
radiative envelope and have a ZAMS chemical composition 
(X,Z) = (0.7,0.02). Since, for main-sequence stars, the in- 
fluence of nonadiabatic effects on the oscillation frequencies 
is known to be small (e.g. Unno et al. 1989), the eigenfre- 
quencies are determined in the isentropic approximation. 

The frequencies of the lower-order <7 + -modes typically 
range from 0.1 T^ n to 10 T^ n , so that the condition that the 
forcing frequencies be in the range of the eigenfrequencies 
can be expressed as 

kn + mQafT^ (10) 

with 0.1</<10. For rotational angular velocities Q <C 
r dyn' the- kih harmonic in the expansion of the tide- 
generating potential may thus give rise to resonances with 
lower-order g + -modes for orbital periods 



Figure 6. As Fig.0 but for TMS stars. 
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Figure 7. The frequency range covered by the eigenfrequencies 
of the second-degree g+-modes of radial orders N = 1, 2, . . . , 25 
in main-sequence stars with a mass between 2 Mq and 20 Mq 
and an evolutionary state between the ZAMS (dashed lines) and 
the TMS (solid lines). The frequencies are expressed in units of 
the inverse of the star's dynamical time scale. 



Poxb ~ — J — Tdyn- (11) 

In order to constrain the forcing frequencies leading to 
significant resonances between dynamic tides and free oscil- 
lation modes, we consider Eq. li lt in the cases where / = 0.5 
and / = 5.0. The former choice represents the bulk of the 
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lower-order g + -modes displayed in Fig. |7| while the latter 
choice represents an upper limit on the frequency range of 
the g + -modes. For presentation purposes, the two values of 
/ are conveniently chosen to differ by a factor of 10. The 
choice is convenient because, for a given stellar model, the 
orbital periods associated with k and / are equal to those 
associated with 10 k and 10 /. The variations of the result- 
ing orbital periods as a function of the stellar mass are dis- 
played by the solid lines in Figs. |H]andQ2] for stars at the 
beginning and at the end of the main sequence, respectively. 
The value of k corresponding to each line is indicated in the 
tables on the right-hand side of the figures. It follows that, 
for ZAMS stars in binaries with orbital periods of the order 
of a few days, the lower-order harmonics in Expansion J3J of 
the tide-generating potential give rise to forcing frequencies 
of the order of ~ 0.5r^, n , while higher-order harmonics give 
rise to forcing frequencies of the order of ~ 5.0t^. For 
stars on the TMS, the corresponding orbital period range 
increases to a few tens of days due to the longer dynamical 
time scales associated with these stars. 

Since the small parameter et appearing in the definition 
of the tide-generating potential is inversely proportional to 
the square of the orbital period, the tidal force decreases 
rapidly with increasing orbital separations. From Eqs. @ 
and 1111 . one derives that the values of et associated with 
the orbital periods represented by the solid lines in Figs. |H1 
and are given by 



et = 



t 

k 2 



1 + 9 



(12) 



For binaries with mass ratio 5 = 1, the parameter et takes 
the value 10" 2 when / = 0.5 and k w 3 - 4 (or / = 5.0 
and k w 35 — 36) and decreases to 10 -3 when / = 0.5 and 
k « 11 - 12 (or / = 5.0 and k « 111 - 112). Smaller mass 
ratios yield smaller values of k. 

As discussed in Section the rapid decrease of the 
Fourier coefficients C2, m ,fe in the expansion of the tide- 
generating potential puts an upper limit on the maximum 
value of k that needs to be considered for a given value 
of the orbital eccentricity e. Vice versa this implies that 
for each value of k there is a minimum orbital eccentricity 
for which the term associated with the Fourier coefficient 
C2,m,fe provides a non-negligible contribution to the expan- 
sion of the tide-generating potential. We determined this 
minimum value for each value of k listed in the tables on 
the right-hand side of Figs. |H| and [§] To this end, we re- 
stricted ourselves to the terms associated with the azimuthal 
number m — — 2 since they generally provide the largest 
contributions. We furthermore adopted two different thresh- 
olds to separate contributing from non-contributing terms: 
|c 2 -2,fc(e)| « 10~ 2 and |c 2 ,_ 2 , fc (e)| w 10" 4 (cfr. the solid 
lines in Fig. 0. The resulting values of e are listed in the 
columns labelled e_2 and e_4, respectively. The table shows, 
for example, that the minimum orbital eccentricity required 
for the absolute value of the Fourier coefficient 02,-2,5 to be 
larger than or equal to 10 -2 , and thus for the associated 
term in the expansion of the tide-generating potential to be 
non- negligible, is e = 0.17. When the absolute value of the 
Fourier coefficient 02,-2,5 is required to be larger than or 
equal to 10 -4 , the condition on the orbital eccentricity is 
relaxed to e > 0.04. 

The high-order harmonics required to reach forcing fre- 



quencies of the order of ~ 5.0 need very high orbital 
eccentricities in order for the associated terms to contribute 
non-trivially to the tide-generating potential. These high or- 
bital eccentricities in combination with short orbital periods 
yield small periastron distances which may cause the star to 
overflow its Roche-lobe when it passes through the perias- 
tron of its relative orbit. The critical orbital periods sepa- 
rating detached from semi-detached systems were derived in 
Section 0] and presented in Figs. |S] and HJ for different combi- 
nations of the binary mass ratio q and the orbital eccentric- 
ity e. For illustration, the orbital periods leading to semi- 
detached systems in the case of the binary mass ratio q = 1 
and the orbital eccentricity e = 0.5 are displayed in Figs. |S| 
and [5] by means of the grey-shaded background. In this ex- 
ample, the critical orbital periods separating detached from 
semi-detached systems increase from ~ 1.5 days for a 2Mq 
ZAMS star to ~ 3 days for a 20 M© ZAMS star, and from 
~ 5 days for a 2 M Q TMS star to ~ 20 days for a 20 M e TMS 
star. The minimum orbital eccentricity e = 0.52 required 
to have a non-negligible contribution to the tide-generating 
potential from the terms associated with fc-values up to 30 
can therefore not be reached by short-period systems with- 
out becoming semi-detached at the periastron. From Figs. 
and|fj] it follows that the same holds true for the other high- 
order harmonics required to reach forcing frequencies of the 
order of ~ 5.0rJ^ n . The situation improves somewhat for 
binaries with mass ratios q smaller than unity, but even for 
q = 0.1 the systems can barely avoid a semi-detached state 
when the stars are located at the periastron of their relative 
orbit. 

Next, we turn our attention to the forcing frequencies of 
the order of 0.5 associated with the lower-order harmon- 
ics in Expansion of the tide-generating potential. These 
harmonics provide non-negligible contributions to the tidal 
potential for low to moderate values of the orbital eccen- 
tricity, so that the occurrence of Roche-lobe overflow at the 
periastron is less of an issue. 

In Section^] we have shown that the order of magnitude 
of any given term in the expansion of the tide-generating po- 
tential is not only determined by the orbital eccentricity e, 
but also by the parameter et which is proportional to the ra- 
tio of the tidal force to the gravity at the star's equator. Since 
the parameter et is the dominant factor in binaries with low 
to moderate orbital eccentricities, we neglect the role of the 
eccentricity in what follows and, somewhat arbitrarily, set 
the limit of what we consider as favourable conditions for the 
resonant excitation of free oscillation modes to et > 10 -3 . In 
the case of a binary consisting of two 5 Mq ZAMS stars, this 
corresponds to an upper limit on the orbital period of ~ 5 
days. A similar upper limit for the occurrence of noticeable 
resonances between dynamic tides and free oscillation modes 
is inferred from the tidally induced radial- velocity variations 
presented by Willems & Aerts (2002) and the apsidal-motion 
calculations performed by Smeyers & Willems (2001). 

The maximum value of k for which the orbital periods 
represented by the solid lines in Figs. [8] and [§] give rise to 
forcing frequencies of the order of 0.5 and for which 
et ^ 10 -3 are obtained from Eq. 1121 as 



15.8 



1 + 9 



1/2 



(13) 



8 B. Willems 



1 2 



10 



Zero— Age Main Sequence 



10 



20 
M,/M Q 




















f = 0.5 






f = 5.0 




— 


k 




e- 4 


k 






- 


I o 




u.zy 


I JU 


>U.o 


> U . o 


- 


1 9 


n aa 


n 9^ 


1 9H 


/U,o 


-> u . o 


- 


1 1 


n a 1 


fl 94. 


1 1 n 


u . o 


X u . o 




1 


0.39 


0.21 


1 00 


>0.8 


0.75 


— 


9 


0.36 


0.18 


90 


>0.8 


0.74 




8 


0.32 


0.1 4 


80 


0.79 


0.72 




7 


0.28 


0.1 1 


70 


0.78 


0.70 




6 


0.23 


0.07 


60 


0.76 


0.67 




5 


0.17 


0.04 


50 


0.73 


0.64 




4 






40 


0.70 


0.59 




3 






30 


0.65 


0.52 




2 






20 


0.56 


0.41 




1 






10 


0.39 


0.21 



















30 



Figure 8. The orbital period range that may give rise to resonances with lower-order g+-modes in 2 — 20 Mq ZAMS stars. The solid lines 
correspond to the orbital periods obtained from Eq. Hill by setting / = 0.5 and k = 1, 2,3, . . . , 13; or / = 5.0 and k = 10, 20, 30, . . . , 130. 
For each value of k, the table on the right-hand side of the figure lists the minimum values of the orbital eccentricity e_2 and e_4 for 
which |c2 _2 k( e )\ ;^10 -2 and |c2 — 2,fc( e )l ^ 10 — 4 ) respectively. The grey-shaded background indicates the orbital periods for which a 
binary with a mass ratio q = 1 and an orbital eccentricity e = 0.5 becomes semi-detached when the stars pass through the periastron of 
their relative orbit. 
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Figure 9. As Fig.|SJ but for TMS stars. 
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The corresponding orbital periods are 

> 1/2 

-forb.max ~ 198.7 Tdyn [ ~1 ! 

\i + q 



(14) 



In the case of a binary with a mass ratio q = 1, the upper 
limits take the values fc max ~ 11 and P rb,max ~ 140rd yn - 
They change by less than a factor of v2 when 1/3 < q < oo, 
but decrease rapidly for binary mass ratios smaller than 1/3. 
The minimum orbital eccentricity required for the terms as- 
sociated with k w 11 to provide a non- negligible contribu- 
tion to the tide-generating potential is e ~ 0.25 — 0.4. The 
longest orbital periods thus still require substantial orbital 
eccentricities in order to be able to excite oscillation modes 
with eigenfrequencies of the order of 0.5 t~[*. 

The results shown in Figs. |H| and [5] can furthermore 
be extrapolated to resonances between dynamic tides and 
free oscillation modes with eigenfrequencies smaller than 
0.5 T^y n . From Eq. 1111 it follows that the orbital periods 
giving rise to such resonances increase with decreasing order 
of magnitude of the eigenfrequencies, so that the condition 
£t > 10 -3 becomes more and more stringent until it can no 
longer be fulfilled for any value of k and q. The smallest 
multiple / of the star's dynamical time scale for which the 
orbital periods resulting from Eq. 11 U yield values of et 
larger than 10~ 3 is given by 



fn 



0.03 



1+q 



1/2 



(15) 



[see Eq. 1)12^ ]. For binaries with a mass ratio q = 1, this 
imposes a lower limit of ~ 0.05 T^ n on the frequency range 
favourable for resonances between dynamic tides and free 
oscillation modes. In the case of a binary with a mass ratio 
q = 0.1, the lower limit increases to ~ 0.1 . The condi- 
tions for the resonant excitation of free oscillation modes g + 
thus rapidly become less favourable with decreasing order of 
magnitude of the involved frequencies. 

So far, we only considered slowly rotating stars for 
which Q <C T^ yn . In more rapidly rotating stars the rota- 
tional angular velocity Q may significantly affect the deter- 
mination of the forcing angular frequencies induced by the 
companion. Resonances with eigenfrequencies of the order 
of / times the inverse of star's dynamical time scale then 
occur for orbital periods 

27rfc 



■ ,. , , -Tdyn- (16) 

/ — mS2r dy n 

For m — 0, the equation reduces to Eq. I|ll[l so that the 
results presented in Figs. [§] and |2l generally remain valid. 
However, since the coefficients C2,o,fc decrease more rapidly 
with increasing values of k than the coefficients C2._2.fc the 
associated terms in the expansion of the tide-generating po- 
tential require a higher orbital eccentricity in order to pro- 
vide a non-negligible contribution. For m = —2, on the other 
hand, the periods represented by the solid lines in Figs. [£] 
and [5] decrease with increasing rotational angular velocities 
fl. Consequently, for both m = and m = —2, the upper 
limit of 5.0 t~[ on the frequency range that is favourable 
for the resonant excitation of free oscillation modes becomes 
even more stringent for rapidly rotating stars than it was for 
more slowly rotating stars. The upper limit on the range of 
orbital periods favourable for the resonant excitation of os- 



cillation modes with eigenfrequencies of the order of 0.5 rT^ 
is not affected by the rotational angular velocity Q, although 
the value of k for which the upper limit is reached increases 
with increasing values of £1. Our qualitative results thus re- 
main valid for rotation rates that are somewhat higher than 
implied by Eq. Q. In the case of very fast rotators, the 
situation becomes significantly more complicated since the 
forcing frequencies may take negative as well as positive val- 
ues. We refer to Willems & Claret (2003) for a more detailed 
account on the behaviour of the forcing frequencies in such 
stars. 

Finally, we note that for a rotating star the frequencies 
and properties of <? + -modes with eigenfrequencies of the or- 
der of the rotational angular velocity Q may be significantly 
altered with respect to those of a non-rotating star. In par- 
ticular, stellar rotation lifts the degeneracy of the eigenfre- 
quencies with respect to the azimuthal number m, so that 
more modes become available which may be resonantly ex- 
cited (e.g. Unno et al. 1989). Stellar rotation furthermore 
induces additional low-frequency modes known as inertial 
modes and r-modes which have frequencies of the order of 
the star's rotational angular velocity Q. The resonant ex- 
citation of these modes may lead to comparably important 
effects as the resonant excitation of free oscillation modes 
g + (Rocca 1982, Savonije & Papaloizou 1997, Papaloizou & 
Savonije 1997). 



6 CONCLUDING REMARKS 

In this paper, we initiated a systematic study aimed at ex- 
ploring the range of stellar and orbital parameters leading 
to favourable conditions for the resonant excitation of g + - 
modes by tides in close binary components. We considered 
slowly rotating stars with masses between 2 and 20 M©, and 
evolutionary stages at the beginning and at the end of the 
main sequence. 

A crucial ingredient in the investigation is the expansion 
of the tide-generating potential in Fourier series in terms of 
the mean Keplerian motion. The expansion induces an in- 
finite number of forcing frequencies in the star which may 
be close to the eigenfrequencies of the star's free modes of 
oscillation. In practice, the rapid decrease of the Fourier co- 
efficients limits the number of forcing frequencies to a finite 
number which increases with increasing values of the orbital 
eccentricity. The number can be derived from the variations 
of the dominant Fourier coefficients shown in Figs. 

For given orbital and rotational periods, the finite 
number of non-trivial terms in the expansion of the tide- 
generating potential associated with a given value of the 
orbital eccentricity puts an upper limit on the forcing fre- 
quencies induced in a star, which in turn puts an upper limit 
on the eigenfrequencies of the oscillation modes that may be 
resonantly excited. The limit is of the order of five times the 
inverse of the star's dynamical time scale and is particularly 
stringent for the /-mode and the lowest-order g + -modes in 
evolved low-mass main-sequence stars (see Fig. 0. It also 
contributes to the difficulties to excite p-modes by tides in 
close binaries (see also Uytterhoeven et al. 2003, in prepa- 
ration) . 

Resonances between dynamic tides and g + -modes with 
frequencies of the order of half the inverse of the star's 



10 B. Willems 



dynamical time scale on the other hand are constrained 
by the rapid decrease of the strength of the tidal force 
with increasing orbital separations. When the limit sepa- 
rating favourable from less favourable conditions is set at 
St ~ 10 , where et is a measure for the ratio of the tidal 
force to the gravity at the star's equator, the highest order 
harmonics in the expansion of the tide-generating potential 
that may give rise to such resonances are given by Eq. I13H . 
This corresponds to an upper limit on the orbital period 
given by Eq. 1141 . For oscillation modes g + with eigenfre- 
quencies much smaller than half the inverse of the star's 
dynamical time scale, the conditions for resonant excitation 
quickly become less favourable. However, these very low fre- 
quencies are only reached by g + -modes of very high radial 
order, so that this behaviour does not impose severe con- 
straints on the parameter space available for the resonant 
excitation of free oscillation modes. 
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